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Abstract 

Several nucleon-nucleon potentials, Paris, Nijmegen, Argonne, and those derived by 
quantum inversion, which describe the NN interaction for T^ab ^ 300 MeV are extended 
in their range of application as NN optical models. Extensions are made in r-space using 
complex separable potentials definable with a wide range of form factor options including 
those of boundary condition models. We use the latest phase shift analyses SPOO (FAOO, 
WIOO) of Arndt et al. from 300 MeV to 3 GeV to determine these extensions. The 
imaginary parts of the optical model interactions account for loss of flux into direct or 
resonant production processes. The optical potential approach is of particular value as 
it permits one to visualize fusion, and subsequent fission, of nucleons when T^ab > 2 
GeV. We do so by calculating the scattering wave functions to specify the energy and 
radial dependences of flux losses and of probability distributions. Furthermore, half-off 
the energy shell t-matrices are presented as they are readily deduced with this approach. 
Such t-matrices are required for studies of few- and many-body nuclear reactions. 



1 Introduction 

A theoretical description of nucleon-nucleon (NN) scattering is a fundamental ingredient for 
the understanding of nuclear structure and scattering of few- and many-body nuclear systems 
im, ^, ^. This is a paradigm of nuclear physics. Of the spectrum, low energy NN scattering 
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traditionally is described in terms of few degrees of freedom of which spin and isospin symmetries 
play the predominant role. At medium energies, production processes and inelasticities become 
important and several elementary systems composed of nucleons and mesons contribute to NN 
scattering. While these nucleons and mesons are emergent structures from QCD, at present 
there is no quantitative description of NN scattering above the inelastic threshold either in 
terms of QCD or of the emergent nucleons and mesons 0]. 

Theoretically undisputed is the need for relativity of which there are two aspects. First is 
the increasing importance of relativistic kinematics as the kinetic energy becomes comparable to 
the rest masses of the scattering particles. Second, particle production is inherently relativistic 
requiring, ultimately, a description in terms of highly non-linear QCD. But that non-linearity 
inhibits a facile QCD explanation of NN scattering. Notwithstanding, there exist hybrid models 
which offset that non-linearity in seeking explanation of the excitation spectra and of the 
scattering of hadrons p, 0. All use heavy valence quarks, with an effective mass typically of 
300 MeV, and massive Goldstone bosons in lieu of massless gluons. They also maintain color 
degrees of freedom. As well there are effective quantum field theories (EQF) that link the 
quark-gluon structure of the standard model to low energy nuclear physics [§, |^. Currently 
these latter approaches are very popular as they may give a foundation and interpretation of 
emergent structures. But like the hybrid models, due to the underlying expansion schemes used 
with EQF, many ad hoc degrees of freedom are involved. 

The experimental NN data and its parameterization in terms of amplitudes and phase 
shifts, are very smooth with energy to 3 GeV |]TU|, |TT], |T^; a feature which supports use of 



the classic approach using a free NN interaction potential. By so doing one uses a minimal 
number of degrees of freedom with again those degrees of freedom being associated with the 
spin and isospin of the total system. Of course, this classic approach sacrifices all reliance 
on substructures. However, the underlying dynamic still reflects its geometric facet by means 
of surfaces and boundary conditions. The success of bag models is direct evidence of the 
crucial role such boundary conditions play with the emergent structures from them being direct 
consequences of QCD confinement. This is further support for our view that an explanation of 
elastic NN scattering need not, if will not, depend explicitly upon QCD details. Only geometric 
attributes of the QCD domain, we contend, are required to explain most data. This view is well 
supported by high energy scattering for which the geometric limits of the S-matrix are reached 
and form factors are defined independent of energy. In the transition region the geometric limits 
are not reached and the factorization schemes [|1^, Q used at higher energies do not apply. 

Of course, in the last decade or so, there have been several theoretical attempts built upon 
boson exchange models to explain NN scattering data to 1 GeV. All such have given but 
qualitative results, often requiring many degrees of freedom even to achieve that qualitative 
agreement and despite explicit inclusion of A and A^* resonances. Optical model studies have 
also been made for medium and modest high energy NN scattering [|15|, |16|, |T^, and they can 
be improved to give a high quality description of scattering at medium energy. 

A high quality fit of on-shell t-matrices by a potential model is very desirable also as it 
facilitates extension into the off-shell domain; properties which are needed in few and many body 
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calculations. In particular, microscopic optical model potentials for elastic nucleon-nucleus 
scattering and bremsstrahlung reactions that give quantitative results, require a careful and 
exact treatment of the off-shell NN t-matrices Furthermore, calculations of such entities 
have shown that it is crucial to have on-shell values of the t-matrices in best possible agreement 
with NN data at all energies. Concomitantly one needs high precision NN data against which 
one can specify NN interactions. 

There are many studies of few and many body problems in the low energy regime T^ab < 300 
MeV and the results have consequences for any model extension above threshold 0. We note 
in this context that significant off-shell differences in t-matrices are known to exist between the 
theoretically well motivated boson exchange models of NN scattering in this regime. It remains 
difficult to attribute with certainty any particular dynamical or kinematical feature with those 
differences. Non-locality, explicit energy dependence and features associated with relativistic 
kinematics are some possibilities. 

In contrast, there is the quantum inverse scattering approach by which any on-shell t-matrix 



can be continued into the off-shell domain |T^. A specific method is the Gel'fand-Levitan- 



Marchenko inversion algorithm for Sturm-Liouville equations. This approach to specify t- 
matrices off-shell is appropriate when the physical S-matrix is unitary and the equation of 
motion is of the Sturm-Liouville type. Such is valid without modification for NN t-matrices 
in the energy regime to 300 MeV. Mathematically, the Gel'fand-Levitan-Marchenko algorithm 
is a method based upon a class of real and regular potentials. In the spirit of inverse scat- 
tering, we generalize that method for non-unitary S-matrices. By that means we generate an 
NN optical model separately for each partial wave. The algorithm we have developed allows 
studies of complex separable potentials in combination with any background potential. The 
background potential can be any of the existing r-space NN potentials. We have not used 



k-space background potentials, such as Bonn-B g], Bonn-CD ||2^ and OSBEP albeit that 
similar analyses can be made with them. 

We limit the reference potential to the well known real r-space potentials from Paris [^], 
Nijmegen (Reid93, Nijmegen-I, Nijmegen-II), Argonne (AV18), and from inversion I^S], 
To them we add channel dependent complex separable potentials with energy dependent 
strengths. For given input data results then, the full potentials are unique. The experimental 
background and motivation for analyses using an optical model is given in Sect. II. A detailed 
description of the theoretical algorithm is given in Sect. III. A discussion of results is given in 
Sect. IV while a summary is given in Sect. V. 



2 Survey of data and motivation for the optical potential 



NN scattering is a long standing problem which has been reviewed often as the database 
developed |TD|, [IT], The low energy data has been analyzed by the VPI/GWU group 

for TLab < 400 MeV, the Nijmegen group ^ with the NN phase shift results PWA93 for 
TLab < 350 MeV, and by Machleidt M giving the Bonn-CD-2000. Of these, the VPI/GWU 
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group has given many solutions for this low energy regime over the years, all of which have been 
listed by Arndt et al. in a very recent publication |jl2| . For their use note that the solution name 



reflects the season and year of their creation although the low energy solutions have names that 
end with 40. Clearly that database has grown rapidly in the last two decades. While the pp 
data now extends to 3 GeV, the np data are limited to 1.3 GeV. Surprisingly, the solutions 
from SM97 to WIOO remain very closely the same and are very stable with regard to new data. 
We have used the solutions SPOO, FAOO and WIOO in our calculations and found results that 
differ but marginally. Thus hereafter in the main we refer solely to the results of calculations 
based upon SPOO. The findings are equally valid for other more recently dated solutions. In 
our practical applications however when new potentials are sought their generation is based 
upon the most current solution 



The VPI/GWU solutions are parameterizations of the elastic channel NN S-matrix. 
They consider 

Si = {I + iK^){l - iK^)-^ (1) 

which inverts to give 

Ki = i{l- Si){l + Si)-^ = ReK^ + ilmK^. (2) 

The real part of this K-matrix is related to a unitary S-matrix (Sq) and therewith phase shifts 
and e are defined by 

„ _ (l+iReKi) _j cos 25 exp 2^(5" i sin 2e exp i{S~ + S~^) 1 , . 

^ ~ [l-iReKi] ~ \ z sin 2e exp i{5- + 5+) cos 2e exp 2i5+ ]' ^ ' 

The absorption parameters and /i relate to the imaginary part of that K-matrix by 

ImK =[ tanV tanp"tanp+cos/i 1 

^ I tanp~ tanp+ cos/i tan^p"^ J ' 

These relations simplify to K = tan 5 + i tan^ p for uncoupled channels. 

In our study, real NN potentials derived from fixed angular momentum inverse scatter- 
ing theory have been used. They have been generated from inversion algorithms predicated 
upon the Gel'fand-Levitan-Marchenko integral equations which physically link to the radial 
Schrodinger equation of a fixed angular momentum. 



ij,{r,k) = k''ij,{r,k), (5) 



where V^(r) is a local and energy independent operator in coordinate space. Substituting 



g(r) = ^^^ + ^V^,(r), and \ = k\ (6) 
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identifies Eq. (||) as a Sturm-Liouville equation 



dx^ 



+ q{x) 



y{x) = \y{x). 



(7) 



There are two equivalent inversion algorithms for the Sturm-Liouville equation, which one 
identifies as the Marchenko and the Gel'fand-Levitan inversion. Both yield principally the same 
solution and numerically they are complementary. The salient features are outlined for the case 
of uncoupled channels. For coupled channels the inversion equations are matrix equations with 
input and translation kernels correspondingly generalized. 

In the Marchenko inversion the experimental information enters via the S'-matrix, S£{k) = 
exp{2i6£{k)) , with which an input kernel is defined in the form of a Fourier-Hankel transform 



1 

2^ J-c 



hjirk) [Si{k) - l]hj{tk)dk, 



(8) 



where h^{x) are Riccati-Hankel functions. This input kernel when used in the Marchenko 
equation, 

/oo 
AfXr,s)F,{s,t)ds = 0, (9) 

specifies the translation kernel A£{r,t). The potential of Eq. (|^) is a boundary condition for 
that translational kernel, 

V,{r) = -2-^A,{r,r). (10) 

The Gel'fand-Levitan inversion requires not the S'-matrix but rather the Jost-function as 
spectral input. The latter is related to the S'-matrix by 



Si{k) 



F,{-k) 
Fi{k) ' 



The Gel'fand-Levitan input kernel then is defined as the Fourier-Bessel transform 



2 

TT Jo 



1 



je(tk)dk, 



where je{x) are Riccati-Bessel functions. The Gel'fand-Levitan integral equation 

K,{r, t) + Ge{r, t) + [ K,{r, s)Ge{s, t)ds = 0, 
Jo 

also defines a translational kernel with boundary condition 



(12) 



(13) 



Ve{r)=2—Ke{r,r). 



(14) 
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The boundary conditions Eq. ( [T0| ) and (|T^ yield identical potentials. 

Determination of the input kernels from data, phase shift functions S{TLab{k)), or K- matrices 
K{TLab{k))i requires an accurate interpolation and extrapolation of that data. In all practical 
applications rational functions are very appropriate. In this work we made a representation 
of data for Tiabik) < 3 GeV where the order of the rational functions p[^^~^'^^l(fc) was 
chosen to be as small as possible. Typically 2 < < 6. An implication is that extrapolations 
of 6{k) from the highest energy (last) data point kmax to infinity do not change sign and 
limfc_,oo S{k) ~ l/fc. We control the rational function fit with weight functions which guarantee 
that those fits will be particularly accurate for some desired interval and less stringent elsewhere. 
For example, the channels ^5*0, ^Pi, ^Po,i, and were weighted with wlow = 1 for 
Tiab < 1-2 GeV and for larger energies, wmgh = 0.05. For the and ^F^ channels, the 
cut between wlow and wmgh was 300 MeV. Consequently the rational functions used in the 
inversion algorithm ensure that the resulting potentials will give the desired values of phase 
shifts from solutions of the Schrodinger equation. Such is evident from the comparisons given 
in Figs. |I] and ^ Therein the fits to the phase shifts to 300 MeV resulting from all three 
models are considered as high quality. Single and coupled channel phase shifts from SM97 
and FAOO solutions for J < 3 are shown together with values found from calculations made 
using three potential models. These model phase shifts were generated with Nijmegen-I and 
Nijmegen-II [E3|, and Argonne AV18 interactions, and with potentials determined using 



Gel'fand-Levitan-Marchenko inversion 17, 25, 26 



On the scale to 3 GeV the OBE model results clearly diverge from data. As with the 
phase shift analyses, OBE potentials have received several critical reviews]^, including 
observations that there are small variations between phase shift analyses and potential model 
results in the subthreshold domain T^ab < 300 MeV A theoretically stable result would 
require many quantities, that need be specified a priori, to be determined from other sources. 
At present that does not seem feasible and all current potentials rely upon fits of many of 
their parameters to the same data. All such fits, however, have been made independently of 
each other and are based upon differing theoretical specifications of the boson exchange model 
dynamics. In Figs. |] and ^ we give a quantitative demonstration of the ensuing differences. 
Therein the Nijmegen phase shift analysis PWA93 |]2^ has been used as reference values for 
various other phase shift solutions and potential predictions for the np ^Sq, and ^Pq channels. 
Such differences are characteristic of variations between finite power series expansions of data 
in a finite interval. A mathematical property of such finite power series expansions within an 
interval is that, while the data in the interval will be well reproduced, continuations beyond 
that interval can radically diverge. Such a property is in evidence in Figs.[l| and 0, and that 
variance is the reason for the caveat often espoused that use of OBEP beyond the fitted energy 
range should be prohibited. Be that as it may, one could expect from a consistent theory that 
such extrapolations, albeit in error, would be the same. Clearly they are not. However of one 
thing we can be sure, the lack of physics with these models lies within the interaction distance 
< 1 fm. The optical model approach we present is an attempt within the frame of potential 
scattering theory to account for and identify such short range properties. 
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Figure 1: Single channel phase shifts for SM97 TLah < 2.5 GeV), FAOO TLah < 3 GeV, and 
reference phase shifts using inversion (In-HH), Nijmegen (Nij-1, Nij-2) and Argonne (AV18) 
potentials. 
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Figure 2: Coupled channel phase shifts for SM97 TLab < 2.5 GeV, FAOO TLab < 3 GeV, and 
reference phase shifts using Nijmegen(Nij-l, Nij-2) and Argonne (AV18) potentials. 
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It is apropos to make a brief remark on the long range character of the NN potential which 
theoretically is identified with OPEP. In the phase shift analysis PWA93 by the Nijmegen 
group and in that of Bonn-CD-2000 of Machleidt, such character is enforced. Indeed that 
precise character re-emerges when either of those phase shift functions are used as input to 
a Gel'fand-Levitan-Marchenko inversion. On the other hand, the VPI/GWU group makes 
no use of OPEP in any of their solutions. Exactly the same quantum inversion of the SM94 
solution does not give OPEP except on average which might be interpreted as signaling the 
importance of nonlocality||26||. 

Despite limitations as discussed above, the OBEP remain the best motivated potential 
models for low energy scattering. They do yield high quality fits to the phase shifts in that 
domain. Such is useful for us in our quest to interpret data with increasing energy. In Fig.^ 
we show an interaction scheme in terms of radial separation that is suitable for low energy 




l<r<10fm OBEP 

NN-out ^ 4^1 m^-* 

Figure 5: Interaction scheme for low energy scattering, < T^ab < 300 MeV 

scattering. This scheme is supported by potentials determined by inversion which reproduce 
the low energy phase shifts used as input to an accuracy \6{exp.) — 6{rat.)\ < 0.25 degrees. 
Such inversion potentials have been made also to follow closely the SPOO real phase shifts to 3 
GeV and these are shown in Fig. ^ They possess a long range Yukawa tail, a medium range 
attraction ~l-2 fm and a strong short range repulsion with an onset at 1 fm. These potentials 
are energy independent so that the long and medium range potential properties diminish in 
importance for kinetic energies above 500 MeV. For projectiles with Tiab > 1-5 GeV essentially 
only the repulsive core of these potentials remains of significance for scattering. Thus inversion 
potentials have also been obtained with the SPOO real phase shifts to 3 GeV using wlow = 0.1 
for Tiab < 1.2 GeV and wnigh = 1 for higher energies, to emphasize the high energy data and 
fix more stringently the short range (< 1 fm) character of the deduced interaction. The short 
range properties of inversion potentials so found are displayed in Fig. |^. Clearly, the ^Sq and 
■^Po.i inversion potentials based upon SPOO real phase shifts which extend to 3 GeV are soft 
core interactions. We neglected in this analyses the np ^Pi channel due to the limited data set 
for Tiab < 1.2 GeV. The higher partial waves are strongly screened by the centripetal barrier 
and so also are not considered here. The core strengths of these ^5*0 and ^Po,i potentials reach 
a shoulder and maximum with a typical value ~ 1 GeV at a radius of 0.3 to 0.4 fm. It is worth 
noting that the shoulder /maximum aspect of the core is a result of flat minima between 1.5 
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Figure 6: Nucleoli- nucleon inversion potentials using SPOO phases 

and 2 GeV in the ^5*0 and ^Po,i SPOO phase shift functions. For higher partial waves, phase 
shift minima lie beyond 3 GeV. As the experimental phase shifts are limited to 3 GeV we have 
confidence in the specified inversion potentials only to about 0.25 fm. The shorter distance 
values reflect only our extrapolation of these phase shifts being limfc^ooi^l^) ~ 

Above 300 MeV reaction channels open and the elastic channel S-matrix no longer is unitary. 
In Fig.|^ we show the gradual increase of the open channels in NN scattering which includes 
resonances as well as single and multiple production thresholds. Only the A (1232) resonance 
has a low energy threshold and a relative small width of 120 MeV. Therefore it is the only 
resonance we expect to be obviously visible in the energy variation of the elastic scattering 
phase shifts. In particular one notices typical variations in the ^D2, ^F^, and ^PF2 channels. 
Otherwise the phase shifts to 3 GeV are very smooth slowly changing functions of energy in 
all channels. Such is a condition for the suitability of a potential model of scattering governed 
by quasi macroscopic geometric entities. In nucleon-nucleus (NA) scattering, entities of that 
ilk are epitomized by the parameters of Woods-Saxon potentials. For the NN case, we have 
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10^ 




used previously |IT6[ a local Gaussian in this same manner, noticing therefrom spin-isospin 
coupling effects more substantial than found with NA scattering. It is also worth noting that 
the absorption in those NN optical potentials for this energy range were not at the geometric 
limit of a fully absorptive disc. Together with the strong spin-isospin coupling, this property 
infers optical potentials that are strongly channel dependent in contrast to the NA case for 
which assumed central and spin-orbit potentials are partial wave independent. 

The plethora of reaction channels that open to 3 GeV, and the requirement of an NN opti- 
cal potential prescription discussed above, mean that it is an interesting but very challenging 
task for any microscopic model to uniquely link QCD substructures to NN scattering phase 
shift functions in analogy to that successful prescription by which NA optical potentials have 
been determined by folding effective NN interactions with nuclear density matrix elements |[18|] . 
Given the success of optical models for NA scattering and as the distorted wave Born approx- 
imation (DWBA) is well known and successful in nuclear, atomic and molecular physics to 
study inelastic reactions, it is intriguing to conjecture use of the DWBA to analyze inelastic 
reaction channels of NN induced reactions. 

In the spirit of visualization of NN scattering shown in Fig. |^, we now include the importance 
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Figure 8: Thresholds for production processes in NN scattering. 



of the reactive and resonant content pictorially in Fig.^. This we consider relevant for 0.3 < 
Tiab < 2 GeV. The upper limit is significant here as we discuss later, but for now it suffices 
that the potential shoulder and maximum seen in Fig.|^ are ~1 GeV. Now we identify some 
specifics in the 0.5 < r < 1 fm range. We conjecture that the two colliding hadrons are retained 
in hadronic states throughout the process. We allow one of the two nucleons to be excited, say 
into a A(1232), while the other remains in the ground state. The excitation may be exchanged 
between the two hadrons as well, and both nucleons may be excited to an intermediate resonant 
state. The production of mesons then can only occur from one or both of the two separate QCD 
entities. The essential feature is that in the energy range, the predominant scattering processes 
are those retaining identifiable hadronic entities. Within an optical potential representation, 
attendant fiux loss equates to a diffuse absorption extending radially to 3 fm and possibly more. 
The bulk of such absorption however lies significantly within 1 fm. 

It requires 2 GeV and more for the two nucleons to surpass the core potential barrier and to 



fuse into a compound system. This is visualized with the scattering sequences shown in Fig. |10 



An objective of our optical model studies is to substantiate this conjecture of fusion and fission 
of resultant compound dibaryonic systems dominating the scattering for this energy regime. 

To describe this developing system for 300 MeV < T^ab < 3 GeV we will use Feshbach 
theory to specify the optical potential[^]. An important feature of that theory is the projection 



14 




r>iofm mSm Coulomb 

l<r<10fm OPEP 




IT 



0.5 < r < 1 fm excitations 




^ P\ meson 

NN-out^l^ y production 

Figure 9: Interaction scheme for medium energy scattering, 0.3 < Tiab < 2 GeV. 



operator formalism with P and Q subspaces, which divide the complete Hilbert space {P+Q) = 
1, into the elastic scattering channel, the P space, and all inelastic and reaction channels which 
are contained in Q space. This theory then assumes a hierarchy of complication in Q space of 
which doorway states are the simplest. Doorway states are characterized to be the only means 
to leave and to return to the elastic channel. Each doorway state in this approach infers a 
complex and separable component in the optical potential with an energy dependent strength. 
If a very large number of doorway states contribute, the effect equates to a local potential 
operator. This was the basis of our previous study 
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2.1 Formal Potential Model 

It is generally accepted that a valid covariant description of NN scattering formally is given by 
the Bethe-Salpeter equation 

M = V + VgM , (15) 

where A4 are invariant amplitudes that are based upon all connected two particle irreducible 
diagrams. This equation serves generally as an ansatz for approximations. Of those, the 
three dimensional reductions are of great use which allow the definition of a potential p^] . 
In particular, the Blankenbecler-Sugar reduction p2i gives an equation very often used for 



applications with NN scattering |3^. This reduction is obtained from Eq. ([I5|) , which in 
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Figure 10: Interaction scheme for high energy scattering, Tiab > 2 GeV 



terms of four-momenta is 

M{q', g; P) = V{q', q; P) + J d^k V(g', k- P) g{k- P) M{k, q- P), (16) 
where the propagator 



Gik;P) 



(2vr) 



(1) 



;|p - ky -M'^ + ie 



(2) 



(17) 



The subscripts refer to nucleon (1) and (2) respectively. In the CM system P = (a/s, 0), which 
is just the total energy E = y/s. In particular, the Blankenbecler-Sugar reduction of the 
propagator Q uses the covariant form 



QBs{k,s) = 
with positive energy projectors 

Aj)(k) 



Sjko) M2Afi)(k)A+)(-k) 
■(27r)3Efe Is-El + ze 



2M 



(19) 
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A^(q', q) = V(q', q) + / ^V(q'. k)- %!^J m(^, q), (20) 



The amplitudes are now expressed with the reduced terms and they satisfy a three-dimensional 
equation 

^ M^ A?,)(k)Af,)(-k ) 
(27r)3 ^^' ' Ek q2-k2 + 2e 
Taking matrix elements with only positive energy spinors, an equation with minimum relativity 
results for the NN t-matrix, namely 



/a k M 

Using the substitutions 



and 



'''' 

V(q'.q)=(|;)'v(q',q)(^y, (23) 
a simplified form of the t-matrix is obtained. It is the familiar Lippmann-Schwinger equation 

(2^^(^''^)^^3k^^(k'^)- (24) 

Of use is an equivalent Lippmann-Schwinger equation for the wave function. Formally, this 
equivalence is proven with the M0ller distortion operator which relates the free wave function 
with the scattered wave and uses the relation between scattering amplitude and potential, 
T(±)$ = V"(](±)$. Finally, we use the equivalence between the Lippmann-Schwinger integral 
equation and the Schrodinger equation so that 

A + ^F(r) - ^/>(r, k) = 0. (25) 



When we identify the potential scale M with the two particle reduced mass 

M = 2/i = 2 26 

mi + m2 

we guarantee consistency with the low energy limit of the Schrodinger equation and use, therein, 
of NN OBE reference potentials. However, a careful and consistent treatment of the M/E 
factors in Eqs. ( p2D and (|2^) is necessary whenever it is important to take relativity into account. 
Minimal relativity enters in the calculation of k"^ by 

2 



s = {mi + m2Y + 2m2TLab=[^Jk^ + ml + ^k^ + mlJ , (27) 

where 

,.2_ miiTlab + '^miTLah) 



For equal masses this reduces to k"^ = s/ 4 — rri^. 
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3 An algorithm for the optical and boundary condition 
models 



We distinguish between three Hamiltonians. They are the reference Hamiltonian Hq, a pro- 
jected Hamiltonian Hpp, and a full optical model Hamiltonian Ti. The first of these, the 
reference Hamiltonian Hq := T + Vq, invokes a given potential Vq for which one can find 
Schrodinger equation reference solutions. The physical outgoing solutions ipo '■= 4'o{'^,^,E) 
of Hq we suppose gives a unitary S-matrix. We assume further that this Hamiltonian is com- 
pletely specified such that evaluation of any quantity, wave function, S-matrix, K-matrix etc. 
is facilitated. The Feshbach projection operator formalism|^T| is used to give the projected 
Hamiltonian, PHqP = Hpp, derived from Hq. We presuppose completeness, P + Q = 1, and 
a finite rank representation of the Q space 

N N 

Q:=J2\<^i><<^i\=Y.\t><i\, (29) 

1=1 i=l 

with the Q space basis functions |$j > interpreted as doorway states. With these doorway states 
we make the link between the QCD and the hadronic sectors; the latter encompassing nucleons, 
mesons and other free particles. Thus we will assume that meson creation/annihilation occurs 
only in the highly nonlinear QCD sector so that Q space wave functions are projections of 
such processes onto hadronic particle coordinates. The third of our Hamiltonians, the full 
optical model Hamiltonian, comprises the reference Hamiltonian Hq and the proper optical 
model potential V. That potential is complex and nonlocal, viz. separable of finite rank, 
H := T + Vo + V(r, r'; Isj, E). Separable potentials are popular representations of NN potentials 
which are designed to serve quite a wide range of purposes P5 |. 



The Schrodinger equation specified with Ti has regular physical solutions := \l'+(r, k, E) 
whose asymptotic boundary conditions we deem to match with the experimental elastic channel 
S-matrix. Specifically, for these experimental S-matrices we have used the continuous solutions 
SPOO from VPI/GWU p8|. The reference potential Vq and separable potential form factors are 
to be specified in detail with any application. 



3.1 Towards a full optical potential model 

To obtain the optical potential on the basis of a given reference potential, we express first the 
solutions of the projected Hamiltonian in terms of the reference Hamiltonian and the a priori 
defined Q space projector. The Schrodinger equation {E — Hq)\iPq >= and its solutions are 
used to express the solutions of {E — Hpp)\ipp >= 0. The latter is equivalent to the Schrodinger 
equation 

_ Hpp - Hqp - HpQ - Hqq)\4jp >= -Hqp\^p > (30) 
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and the Lippmann-Schwinger equation 

l^p > = > -—-L—.HQp\i;p > 

= \iI;o>~Y.G^\J><J\Hqp\^p>- (31) 

3 

These equations are still very general and do not depend upon a specific representation. 
However, in the following we assume a partial wave expansion in terms of spherical harmonics, 
spin and isospin state vectors and radial functions. The following equations are identified as 
radial equations with the set of quantum numbers suppressed. 

Projector orthogonality PQ = QP = implies that 

=< i\^p >=< 2|^o > - < t\G+HQp\^p >, (32) 

and thus 

N 

< j\HQp\i;p >= ^{< $|G+|<I> >}-/ < zlV^o > • (33) 

i 

The solutions of Eq. ( |3lD can be written in terms of {ipo > as 

N N 

\^P > = 1^0 > - EG'+K > {< -flG+l^ >}^'<Mo >= iV^o > -EG'+A,,|^o >, (34) 

ij ij 

wherein one can identify a separable potential 

\i> {<^\G+\^>}r^' <j\ = \i>Xi,<j\=:Ai,ir,r') (35) 

Note then that definition of Q space gives a specification of the separable strengths Xij{lsj, E) 
that is unique. The resultant Eq. ( p4|) has the form of a first order Born approximation but in 
fact it is an exact result. 

To proceed, we initially abandon the exactitude of Eq. (0) and require the strength matrix, 

A., = {< <1.|G+|$ >},r.\ (36) 

to be constrained asymptotically by the experimental S-matrix of the full Hamiltonian 
Schrodinger equation, i.e. asymptotically we induce iV'p >= >■ This implies that complex 
optical model strengths Xij emerge as a result of matching to Riccati-Hankel functions and non 
unitary S-matrices with 

>= \^Pp >^ 1 [-h'{rk) + h^{rk)S{k)] . (37) 
The strengths A^ then can be simply determined from the linear system of equations 

^h+{Rk) [Sik) - So{k)] = ^ G+|z > A,,- < j#o+ >. (38) 
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To reinforce a Lippmann-Schwinger equation, with the experimental S-matrix as boundary 



condition or equivalently with strengths from Eq. (|3J 
potential Eq. ( |55| ) is made. This is achieved with 

1 



V(r, r') 



A- 



a transformation of the separable 



(39) 



;i - G+A)' 

which contains the separable potentials as defined with Eq. (^5]) but whose strengths now are 
solutions of Eq. (^). As the transformation Eq. ( pOj) contains integration of orthonormal func- 
tions, only strengths are altered. Using this optical model in the full Hamiltonian, physical 
solutions are obtained with reference solutions > and Greens function of the reference 
Hamiltonian Hq by means of the Lippmann-Schwinger equation 



(40) 



3.2 Technical details 

The partial wave radial wave functions of the reference potential satisfy equations 



<(r,fc) 



+ 1) 2/i 



K(r) 



K(0 



1 + 2H(r) 



2Kfrl 



1 + 2VJr) 



Ua{r,k), (41) 



wherein we identify the complete set of quantum numbers by the subscript a. These equations 
we have solved numerically for uncoupled and coupled channels using a Numerov method. 
The potentials Va,Vb,V^ are dependent on the quantum numbers {£,s,j) and are taken from 
the Paris, Nijmegen, Argonne, and inversion r-space potentials as one wishes. The Paris and 
Nijmegen-I are momentum dependent potentials with Vb ^ 0, while the Nijmegen-II, Reid93, 
AV18, and inversion potentials all have Vf, = 0. The physical solutions are matched asymptot- 
ically to Riccati-Hankel functions 



2i 



and normalized by 



-h-{rk) + hl{Tk)Sl{k) 
u+{r, k) 



;i + 2H(r) 
The irregular outgoing wave Jost solutions 

J+{r,k)^hl{rk) 



(42) 
(43) 

(44) 



are calculated in the same way as the physical ones and they define the reference potential 
Green functions by 



Gt{ry,k) 



(2/xM2)i^+(r,A;)J-r(r',fc), r < r' 
(2/iM^)ij-+(r,A;)^+^(r',A;), r > r'. 



(45) 
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where the transpose matrix is signaled by the superscript T. At the asymptotic matching radius 
R 

'^t{R,k)=ijt{R,k)+ Gt{R,n,k)<^^{n)dnK{k) $.(r2)V^^(r2, A;)rfr2, (46) 

JQ Jo 

and taking the difference between the reference and full S-matrix, this reduces to 
^UR,k)~^UR,k) = l,h+{Rk)[SM-Sl{k) 

X) POO 

GtiR,n,k)^^{n)dnX^{k) $„(r2)^+(r2, A;)dr2. (47) 

Jo 

A linear expression for the potential strength Xa{k) results. The strengths are transformed by 
Eq. ( p9D to give final separable potential strengths 

/ roo roo \ —1 

a^{k) =[1- K{k) $«(ri)G+(ri, r^, A;)<D,(r2) dr^^ K{k). (48) 



These strengths (Ja{k) define the proper optical model of Eq. (^Oj), given for the more general 
coupled channel and rank < 3 separable potentials, to be 



V(r,r') = 1$, > a, < $„| = p$„W„C (49) 



where 



$ ._ { *]-i/2(0 $|-i/2(r) ^U/,{r) 1 

" ■ I '^]+i/2(0 *?+i/2(0 j ' ^ ^ 

and the symmetric strength matrices are 

Wo,{k) := ReWij + ilmWi^j = {h^/2fx) aa{k), for i, j" = 1 ■ ■ ■ 6. (51) 

For single channel and rank one potentials, this representation is obviously reduced. 

There are several options one may consider for the separable potential form factors $Q,(r). 
First, any finite rank potential may be chosen with the strengths Xa{k) determined from data at 
several energies around a mean energy. In practice, using a rank > 1 option has been successful 
for single channels but inherent lack of energy dependence for coupled channels strongly favors 
restricting potentials to be of rank one. Next is the choice of radial form factors. As rank one 
potential form factors we have used a.) normalized harmonic oscillator radial wave functions 

= ^i{r,hu!), b.) normalized Gaussian functions $q = Ai'oexp— (r — ro)^/aQ with vq and 
Qq being parameters, c.) a normalized edge function $a(ro) = l/2/i, $Q,(ro ± h) = l/Ah and 
$q('^, «) = otherwise, and d.) a boundary condition model realized by $Q,(ro) = l/h and 
^a(^) = otherwise. The last option is suitable for a sudden transition from the hadronic 
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domain into the QCD domain and back. Of course these are but examples and others may be 
inspired by more exphcit considerations of QCD. 

Solutions of the full problem Lippmann-Schwinger equation, Eq. (^Op, have been generated 
with reference potential solutions and Green functions as per Eq. (^5]) and with separable po- 
tentials whose strengths (Ta{k) are given by Eq. (48). These solutions are readily found from 
systems of linear equations, for single and coupled channels, using a trapezoidal integration 
rule for Eq. (^OD recast as 



oo noo 



^a{r)=^Z{r)+ Gt{r,ri)^a{ri)dnWa{k) ^a{r2)K{r2) dr2 . (52) 



However there is a faster method by which solutions Eq. ([46D as well as half off-shell wave 
function solutions and t-matrices can be found. This we consider next. 



3.3 Evaluation of the half off-shell t-matrix 

The calculation of NA optical potentials for Tiab < 3 GeV requires half off-shell NN t-matrices 
for on-shell k- values k < 6 fm~^, and a correspondingly large range of off-shell values. In 
principle, in applications the off-shell k-values (later identified with q) are needed in integrals 
from — > oo but a reasonable upper limit is g = 2k. A fast and stable method of evaluation 
of such half off-shell t-matrices, when r-space potentials are chosen, is an extension of the 
Schrodinger equation as an inhomogeneous differential equation. The method follows that of 
Van Leeuwen and Reiner |p6|| . 

The most general potentials in our study contain momentum dependent, local and separable 
complex potentials for both single and coupled channels. In particular for the results shown, a 
rank one separable potential with a radial harmonic oscillator form factor ^i[r,huj) has been 
used. The Schrodinger equation then can be cast as 



M{r)—- M{r)^^ - K(r) + K'W + 2Vl{r)- + k' 



Mr,k,q) 



$^(r)AKr) / (^e{x)Mx. q) dx + (^ - q')]i{rq), (53) 







with M.{r) = (1+2 Vb{r)). The regular solutions of which not only must vanish at the origin 
but also asymptotically must match to 

lim #'°)(r, k, q)Afi = ji{rq) + /if '°)(rfc)f '"^(A;^, fc, q) (54) 

to determine the half off-shell t-matrix T^'^''^\k'^,k,q) and the normalization A/"^. Spherical 
Riccati functions are symbolized by jiix), hf{x) and h^x) = ni{x). In the following we 
suppress the channel subscript I as the expressions hold for single and coupled channels. The 
on-shell t-matrix gives the S-matrix by the relation 

S{k) = l + 2iT'^+\k^,k,k). (55) 
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To solve for coupled channels ^SDi, ^PF2, etc., two linear independent regular solutions are 
calculated and Eqs. (148|) , (^) and (|5^ are to be understood as 2 x 2 matrix equations. 

The regular solutions are readily found numerically as follows. First, a regular solution of 
the reference potential Schrodinger equation 

/^(r, k) - V{r)fo{r, k) + A;Vo(r, k) = (56) 

is calculated. Therein V{r) implies all the local potential terms including the centripetal barrier. 
Then a regular solution of the full potential Schrodinger equation with the reference potential 
V{r) and separable potential, 

/{'(r, k) - V(r)/i(r, k) + k'^r, k) = ^r)a{k) < $|/i >, (57) 

is obtained from a particular solution of 

g'iir, k) - V{r)gi{r, k) + ^giir, k) = ^{r)a{k)T, (58) 

where we use =< <l>|/o >, and fi{r,k) = fo{r, k)A + gi{r, k). The factor (matrix) A is 
determined from 

< $1/1 >=< $1/0 >A+< ^\gi > (59) 

and 

A = l- J^-^ <(^\gi> . (60) 

Finally the regular solution /i(r, k) can be multiplied with any complex number (matrix) to be 
a general regular solution of Eq. ([571). 

The half off-shell t-matrix is related to the regular half off-shell wave function ?/'(r, fc,g), 
which satisfies the inhomogeneous Schrodinger equation 



i){r, k, q) = $(r)(T(A;) < > +{k^ - q^)j{rq). (61) 



Asymptotically this wave function is 

^lj^^''^\r,k,q) ^ j(rq) + h^^'''\rk)^T^^^''\k^,k,q). (62) 

k 

A general regular solution of Eq. ( |57D and a particular regular inhomogeneous solution of 
Eq. ( p5D then is needed to satisfy the boundary conditions given in Eq. (^). A particular 
solution of Eq. (|6T|) is obtained in two steps. First, with 



^ =< $1(^2 >=< $1/1 >=< $1/0 >, (63) 
a particular solution is given by 

f2{r,k,q) = J\{r,k)B + g2{r,k,q), (64) 
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where B is determined from 

S = 1 - J^-^ < $1^2 > • (65) 
The off-shell wave function matches asymptotically as 

^(±'°)(r,fc,g) = fi{r,k)Af + f2{r,k,q) ^ j{rq) + h'^^''^\rk)lT{e,k,q). (66) 

The normalization J\f and t-matrix T^'^'^\k'^, k, q) are readily evaluated from the quasi Wron- 
skians 



^T'~^''\k',k,q) = W-'[j,h^^'%W[j,h]Ar + W[j,f2 



(67) 



where we define 



W[a, b] := ^ ; -bn - an- 



at two asymptotic radial points r„_i and r„ = r„_i + /i. The quantities a and 6 can be either 
scalars or matrices. 

It is very convenient to use the Numerov algorithm to solve Eqs. (0), (|58D, and (|6TD. But 
to do so for Eq. ( |53D requires equations without first derivative terms. The above can be made 
so by use of a factorization 



ipir,k,q) = f{r,k,q)V{r), with V{r) = 



1 



1 + 2V,(r] 



(69) 



The resulting equation for f{r,k, q) is 



f"{r,k,q) 



■ + l)/r^ - V{r)k^V{r) + V{r)Va{r)V{r) + lv{r)Vl{r)V( 



fir,k,q) 



+^{r)V{r)(r{k) < V<^\f > +{k^ - q^)ji{rq)V{r) 



(70) 



3.4 Numerov Algorithm 

The solution of radial Schrodinger equations is certainly not new and generally deserves no 
mention. Here, we dwell upon the details since we found the specified elements to have a 
normal form of related problems in other fields of physics and engineering which were tested 
with parallel computing facilities. The Numerov algorithm has been widely used for single and 
coupled channels Schrodinger equations since it gives sufficient numerical accuracy with mini- 



mal operations ||3^ . The standard form of linear homogeneous or inhomogeneous Schrodinger 
equations which we have to solve is 



f:'ir) = Vi,ir)f,{r) + W,ir), 



(71) 
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where Wi{r) = for homogeneous equations. The terms Vij{r) and Wi{r) are easily identified 
in Eq. (|70|). For single channels the algorithm is 

fn+l = "^fn - fn-1 + — (^n+l + + Mn-l) (72) 

or 

= (^2 + i^K^ /n - - ^K-l) /n-1 

+ ^{Wn+l + WWn + Wn-l) • (73) 

These expressions generalize for coupled channels using standard vector and matrix algebra. A 
significant reduction of operations is found by using the substitution 

= (l - fn (74) 

in Eq. (j7|). It gives 

^n+l = 2.^n — C,n-1 + ^n, (75) 

and the inhomogeneous equation 

^n+l = 2e„ - ^n-l +^n + ^ (W^+l + lOW^ + W^-l) , (76) 

with 

= T^in. (77) 
J- 12 

Back-transformations from fi use either of the two possibilities 

/i = 6 + Y^^i, or fi = — . (78) 

4 Properties and discussion of the optical model. 

A range of optical potentials have been generated using the algorithm developed above. As 
reference potentials, the Paris, Nijmegen, Argonne, and inversion potentials have been used. For 
the separable potential form factors, normalized harmonic oscillator functions (HO), ^i{r, hu), 
with 200 < hu < 900 MeV have been used. The same hu is used for all partial waves however. 
For single channels all quantum sets with J < 7 were included while those for J < 6 were 
used with the coupled channels. A superposition of several HO functions with radial quantum 
numbers n = 1, 2, 3 was allowed and with data in intervals T^ab ± 25, TLab =t 50, and Tji^at ± 100 
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MeV, optimal solutions of Eq. ( |i7| ) found using a least square linear equation routine from a 
scientific subroutine library (NAG-library). This procedure was used to determine a single huj 
for energies within 0.5 < T^ab < 2 GeV with an overall low x^- That optimal value is hu = 450 
MeV. For higher energies 2 < T^ab < 3 GeV and low partial waves, this optimal oscillator has 
bound states embedded in the continuum, but as such they are of no concern in this analysis 
and so we used hu = 450 MeV for all energies 0.3 < Tiab < 3 GeV. With rank one separable 
potentials, the HO functions (radial quantum number n = 1) are 



$£(r, /ici;) ~ r^^^ exp — (r/ro)^, with tq 



2r 



(79) 



and with huj = 450 MeV, vq = 0.61 fm. Then with the separable form fixed, it is trivial to solve 
Eq. ( ^71) with S-matrix data taken at each energy. In Figs. and ^ we show the full potential 
model phase shifts that result on solving scattering from the deduced optical potentials. They 
are identical with the SPOO solution. 

The strengths (Ja{k) of Eq. (|^) were determined independently for each given reference 
potential and the optical potential values of Eq. (^) are shown in Figs. 13 and H. The lettering 
in the small sub-figures identify the channel. The curves give the results obtained when the 
Paris (Pa), Nijmegen-I (Nl), Nijmegen-II (N2), Argonne AV18 (Av) and single channel inversion 
potentials (In) were used as reference potentials. These optical model strengths display two 
most important features. The first is that they are not insignificant. The reference potentials by 
themselves fail to account for the phase shifts 5 and 5^, e. The second feature of importance is 
the loss of unitarity of the S-matrices accounted for by p and p^, /i. The two features are weakly 
coupled by the optical and reference potentials respectively. Below threshold however, a purely 
real optical potential and very small strengths reflect the agreement of the reference potential 
phase shifts with SPOO. The imaginary potentials show a smooth energy dependence starting 
at threshold Tiab = 280 MeV and, by having negative values, account for flux loss. Notice also 
that the results using inversion reference potentials (In) in the channels ^5*0, ^Pq, ^Pi, ^Pi, 
and have small values for the optical potential real strengths. Thus those real potentials 
need hardly any modifications at short distances and this supports the conjecture of the soft 
core potential discussed in Sect. II regarding Fig.|^. All reference potentials are most uncertain 
in the ^D2, ^F^ and ^PF2 channels. This is well known as the region 300 MeV to 1 GeV is 
dominated by the A(1232) resonance while many A^* and higher spin resonances shape the 
region 1 to 2 GeV. Indeed the obvious energy dependences seen in the and ^F^ channels 
are signatures of the strong coupling to the A(1232) resonance between T^ab = 500 and 750 
MeV. The coupled channel results shown in Fig. |14| follow closely the conclusions drawn for the 
single channel results. Thus only the ^D2 and ^F^ channels show energy dependences in the real 
phase shifts 6 and absorptions p that require particular attention and an explicit treatment of 
resonance coupling. The ^PF2 coupled channels show some similar A(1232) resonance coupling 
around 600 MeV. All the other channels support an energy independent local reference potential 
which can be generated by Gel'fand-Levitan-Marchenko inversion using the real phase shift 
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Figure 11: SPOO phase shifts for np single channels. 
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Figure 12: SPOO phase shifts for np coupled channels. 
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Figure 13: np single channel separable potential strengths, using inversion (In), Paris 
Nijmegen (Nl, N2) and Argonne AV18 (Av) as reference potentials with hw = 450 MeV. 
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Figure 14: np coupled channel separable potential strengths, using Paris (Pa), Nijmegen 
N2) and Argonne AV18 (Av) as reference potentials with hu; — 450 MeV. 
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data. Also, as the optical potential strengths vary smoothly with energy for these channels, 
use of a complex but local very smoothly energy dependent complex potential with Gaussian 
or Yukawa form factors is suggested |T^. It may be that within QCD hybrid models such a 
local background optical potential can be formulated microscopically and be linked with the high 
energy diffraction and Regge models of elastic scattering [0, |13|, |14[ 



In addition to those optical model potentials found by using fiu = 450 MeV, calculations 
where also made using hu = 750 and 900 MeV. This increase in hu reduces the range tq 
from 0.61 0.47 and 0.43 fm respectively. The primary purpose of those calculations 
was a search of the effective radial domain in which the reference potentials all differ most 
markedly. A shorter range of the form factor leads to increased values of the optical 
potential strengths and thus the shortcomings of the reference potentials are magnified. To 
interpret this magnification, one must bear in mind the boundary condition on wave functions 
to be zero at the origin and the influence of a potential in the short range region, between 
< r < 0.8 fm for the results found using huj = 450 MeV and < r < 0.5 fm in the case of 
huj = 750 MeV. Of note in these calculations is that only on using the inversion potentials as 
reference do the real optical model strengths remain small. Given that the inversion potentials 
were designed by themselves to give the SPOO real phase shifts up to Tiab = 3 GeV as derived 
from the real parts of the K- matrix Eq. (j^), that aspect lends further support for a decoupling 
of the real and imaginary parts of the optical model potentials in calculations. Interference 
effects are small with the implication that the real and imaginary parts might be independently 
assessed. Such is not so evident when the OBEP are used as the reference potentials and the 
particular poor extrapolations one finds on using the Paris and Nijmegen-I that have explicit 
momentum dependences are most noticeable. 

The and channel results are exceptional. Even with the inversion potentials as 
reference, the modulating optical potentials have comparable real and imaginary parts. Such 
reflect the means by which the optical model accounts for specific strong resonance effects. 

The changes wrought in complex potential correction strengths when any OBEP is used as 
reference and when the hu for the defining optical potential correction form factors is enlarged 
to 750 and 900 MeV respectively, further stresses that any are poor choices as reference as one 
forces their propriety to even shorter ranges. In sum such have scant credibility in the range 
0.5 < r < 1 fm. 

Our studies support the conceptualization of the formation and fusion of two nucleons, more 
generally of two elementary particles like vrA^, tttt, etc., into a combined object |T^. Such is 



correlated with selective enhancements of probability density and with loss of flux from the 
elastic scattering channel. The probability density of the full problem is 

p„(r, k) = ^ Trace ^^^^(r, /c)^„(r, k) (80) 
and the flux loss function, which results from the continuity equation dtpaij) + (V ■ j)^ = 
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and the time dependent Schrodinger equation, is 
i 1 



(V ■ j)„ = -—Trace / {^l(r, A;)V„(r, ri)^«(ri, k) - ^l{n, A:)Vi(ri, r)^«(r, k)}dn. 



For several low partial waves, in Fig. |T5| we show probabilities as defined by Eq. ( pO]) and 
fiux loss via Eq. (|8TD. In this figure the SPOO phase shift functions diTiab) and piTLab) are 
given as well for each channel and they are compared with the scatter of single energy solutions 
of SPOO. The inversion potential phase shifts are given as well. We show these single energy 
solutions to acknowledge their scatter about the smooth SPOO solutions. Those sharp variations 
have been considered |jll]| as evidence of narrow dibaryon resonances. Should they be so, we 



contend that our potential model and associated viewpoint of fusion is still appropriate on 
geometric grounds. Such dibaryonic resonance effects require a detailed QCD description of 
their structure and decay. 

The contour plots give the probability distributions and the zonal fiux losses for < r < 2 
fm and T^ab < 3 GeV. From these contour plots we envisage a smooth development with energy 
for scattering in all channels with possible exception of the ^Pi, and ^F^ channels. Of those, 
the ^Pi channel is bound by data only to 1.2 GeV, above this value the SPOO phase shift function 
is conjecture. Nevertheless we have used the solution to demonstrate what implication such a 
drastic variation of 5{TLab) for 1 < T^ab < 2 GeV causes in the probability distribution leaving 
the fiux loss essentially invariant. 

The ^5*0 and ^Pq results are given in Fig.|l^. They have very similar characteristics. The 
SPOO continuous energy solutions have phase shifts whose real parts have a minimum at about 
1.6 GeV. The probability and fiux loss plots show characteristic strongly distorted structures 
with the short distance 0.25 < r < 0.6 fm attributes indicative of a large width (P < 1 GeV) 
resonance with strong absorption. The ^Pi results given in the bottom half of Fig. |T3| are 
interpreted similarly. The ^Pi results shown in this figure have more variation as the resonance 
impact in the SPOO solution is refiected in the fiux loss plot in particular. The and ^F^ 



channel results are given in Fig. |T5|. Concomitant with the structured SPOO phase shift functions 
the probability plots indicate a change from the characteristic smoothness of the other channels 
with notable features for 400 < TLab < 900 MeV. A very long ranged probability peak with 
strong distortions and significant absorption extending beyond 1 fm is evident. 

The details shown are not independent of the chosen geometry of the optical potential but 
the patterns are quite stable with variations of the HO energy. These results support our 
pictorial conjectures of reaction schemes given in Figs. ^ ^ and [1^. The energy dependence 
indicates that in the energy regime 300 MeV to 1 GeV the concept that one or the other of the 
colliding hadrons at most is excited to form the A-resonance while the two hadrons remain as 
disparate entities. At higher energies, and for smaller radii, the strong absorption is consistent 
with a fusion of the colliding particles. 

The Kowalski-Noyes f-ratios of the half off-shell t-matrices 

Ti^''Knab{k).Kq) 

Tt''\TLab{k),k,k) 
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Figure 15: Block matrices containing 6{TLab) and piTiab) of SPOO, and energy versus radius 
the probability and flux loss. 
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are useful quantities as they stress the potential differences in momentum space. For a purely 
real potential the Kowalski-Noyes f-ratio is real but this is no longer the case for complex 
potentials. Nevertheless, the f-ratios are always independent of the boundary conditions used 
in Eq. (^) to determine T^^'*^)(fc^, k, q). We show in Fig. ^ a contour plot of the and 
channels for T^ab from 300 MeV to 1 GeV and off-shell momenta q = kojf from to 7 fm~^. 
The Nijmegen-II and inversion reference potentials are used with these calculations. 

The and ^F^ channels were selected specifically as they are noticeably influenced by 
the A (1232) resonance. They also have the most drastic variations of optical potentials with 
the choice of reference potential. The results support our expectations, associated with strong 
energy dependences and/or large differences of experimental and reference potential phase 
shifts, which led to a scattering scheme shown in Fig. ^. It is not difficult to foresee great 
problems in microscopic analyses which attempt to describe the interferences between back- 
ground and resonance scattering, and which aim for a unique high quality result. 

For energies above 1 GeV no obvious resonance effect can be identified with elastic scat- 
tering phase shifts. However, this smoothness does not imply that the off-shell t-matrices are 
independent of the choice of optical potential parameterization. In Fig. |l^ we show in the com- 
plex plane several Kowalski-Noyes f-ratios for the ^Pq channel. In three cases we used HO form 
factors with hu = 450, 750 and 900 MeV, and in one case we used a normalized edge function 
{tq = 0.45 fm, h = 0.015 fm) of Sect. III.B as boundary condition. Quite similar results were 
found for the other channels and the off-shell differences between these results are significant. 
But the influences of such large and obvious off-shell differences disappear when those off-shell 



t-matrices are used in few- and many-body calculations |T8[. It is generally argued that only 
near on-shell values enter in few- and many-body calculations and symmetric sampling around 
the on-shell point implies that any effects of such differences are annulled. Thus we do not 
expect medium energy few- and many-body calculations to be more revealing than were the 
results of calculations at low energy. We consider it not opportune to seek or nominate a 
preference for any of the off-shell t-matrices or particular form factors. 



5 Summary and conclusions 

Diverse nucleon-nucleon r-space potentials, that yield quality fits to NN scattering phase shifts 
for energies below 300 MeV, have been extended to be NN optical potentials from which the 
SPOO phase shift functions to 3 GeV are matched. Complex short range separable potentials, 
addressed as the optical model potential and distinguished from the real reference potentials, 
bridge the gap between the experimental and reference potential phase shifts. By extending 
boson exchange motivated NN potential models to be optical models we invoke a new reaction 
scheme. At medium energy, 300 MeV to 1 GeV, this approach identifies intrinsic excitation of 
isolated nucleons without their fusion. At higher energies, and in particular for energies Tiab > 2 
GeV, the two nucleons can fuse into a compound system, from which meson production and 
other reactions eventuate, as can condensation back into the elastic channel. This view is based 



34 



Re(tjt ) In-HHb (np) 

' off on' 2 




""(Won) '"-HH'd^("P) 




0.4 0.6 0.8 




'■"(Won) '^''-2 °2("P) 




Re (t „/t ) In-HH (np) |m (t „/t ) In-HH (np) 

* off on' 3 ' off on' 3 




0.4 0.6 0.8 1 0.4 0.6 0.8 1 




Figure 16: Kowalski-Noyes f-ratios, real part in left column and imaginary part in right column, 
for the and ^F^ channels calculated with inversion-HH and Nijmegen-II reference potentials 
and optical potentials using HO, huj = 450 MeV, separable form factor. 
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Figure 17: Kowalski-Noyes f-ratios for the '^Pq channel calculated with Nijmegen-II as reference 
potential and optical potentials using HO fiuj = 450, 750, 900 MeV, and a normalized edge 
function (BCM) separable form factor. The dots on the curves fall onto integer koff values and 
kon is that for TLab — 2 GeV. 
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upon the character of the ^5*0 and ^Po,i partial wave phase shifts. Notably it is the minimum in 
the real phase shifts of these channels which transform into soft core potentials. The reaction 
volume of the fused system fits well within a sphere with radius 1 fm and the medium and 
long range boson exchange contributions are small corrections at best. While data at even 
higher energies may indicate a similar reaction scheme with the higher partial waves, it must 
be borne in mind that the centripetal barrier screens that scattering so reducing markedly the 
probability of fusion. 

In the 300 MeV to 1 GeV regime, the A- resonance dominates ^D2, ^F-^ and '^PF2 partial 
waves and all reference potentials require large and strongly energy dependent contributions 
from the optical potential. Our results complement the view that this resonance must be treated 
explicitly. In our case, the A generates a doorway state to pion production and should be 
treated as such within the NN potential model generalization. The separable optical potential 
was chosen to accommodate doorway state formation and decay within a small energy region. 

The OBE reference potentials presently available either give results too far from reality to 
qualify as background phase shifts or use the A-resonance in a way that prohibits separation 
from the background. However, by dint of their construction, inversion algorithms will help 
resolve these issues. The approach is such that one may start with any desired phase shift func- 
tion as input. Of these any real part may be taken as the reference potential phase shifts, whose 
use as input to Gel'fand-Levitan-Marchenko inversion give the reference potentials themselves. 
Therewith, the inversion algorithm we have developed herein can then be used to determine 
the remaining parts of the full NN optical potential. This algorithm facilitates specification not 
only of complex separable potentials, appropriate for specific doorway state effects, but also of 
local complex potentials which encompass smooth energy dependent processes that contribute 
to medium to high energy NN data. The geometric attributes of the optical model, in partic- 
ular the inherent soft core nature of potentials, thus have been determined solely from data. 
Detailed interpretation of these emergent results, of course, must eventuate from QCD inspired 
models. 
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